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We investigate the effect of elastic forward scattering on the ARPES spectrum of the cuprate 
superconductors. In the normal state, small angle scattering from out-of-plane impurities is thought 
to broaden the ARPES spectral response with minimal effect on the resistivity or the supercon- 
ducting transition temperature T^. Here we explore how such forward scattering affects the ARPES 
spectrum in the d-wave superconducting state. Away from the nodal direction, the one-electron 
impurity scattering rate is found to be suppressed as u approaches the gap edge by a cancellation 
between normal and anomalous scattering processes, leading to a square-root-like feature in the 
spectral weight as tJ approaches — Ak from below. For momenta away from the Fermi surface, our 
analysis suggests that a dirty optimally or overdoped system will still display a sharp but nondis- 
persive peak which could be confused with a quasiparticle spectral feature. Only in cleaner samples 
should the true dispersing quasiparticle peak become visible. At the nodal point on the Fermi 
surface, the contribution of the anomalous scattering vanishes and the spectral weight exhibits a 
Lorentzian quasiparticle peak in both energy and momentum. Our analysis, including a treatment 
of unitary scatterers and inelastic spin fluctuation scattering, suggests explanations for the some- 
times mysterious lineshapes and temperature dependences of the peak structures observed in the 
Bi2Sr2CaCu2 08 system. 

PACS numbers: 74.72.-h,74.25.Jb, 74.20.Fg 



I. INTRODUCTION 

Since the earliest ARPES studies of optimally doped 
Bi2Sr2CaCu208 (BSCCO), spectral features near the 
Fermi level have been reported whose width suggested 
the existence of a significant elastic scattering contri- 
bution which varied along the Fermi surface, taking its 
smallest value along the (tt, tt) diagonal and largest near 



the (tt, 0) regioni 



1.2. 3.4.5 



The shapes of the measured en- 



ergy dispersion curves (EDCs) are quite unusual, sug- 
gesting that the 1-electron self energy is also strongly 
energy dependent. Abrahams and Varmai& attempted to 
account for both features by assuming that the self energy 
was a sum of an energy-independent part arising from 
small-angle (forward) scattering, and a second, momen- 
tum independent term modelled with marginal Fermi liq- 
uid theory. They noted that the elastic forward scatter- 
ing was most probably associated with impurities which 
were located away from the Cu02 planes and reflected 
the vf{}^)~^ variation of the momentum- resolved density 
of states as k moved around the Fermi surface. They 
further remarked that, despite the large scattering rates 
of order lOOmeV deduced from fits to ARPES spectra, 
the forward scattering nature of the disorder would be 
consistent with such scattering having a negligible effect 
on the resistivity of the optimally doped cupratesi. In 
addition, it was shown that this type of disorder would 
have a small effect on T^. 

There are two obvious difficulties with this scenario. 
The first is that the spectral peak measured by ARPES 
near the (tt, 0) point is known to sharpen dramatically 
when one goes below Tc, a phenomenon interpreted as 
the formation of a coherent quasiparticle in the super- 
conducting state. This sharpening has normally been at- 



tributed to the well-known collapse of the inelastic scat- 
tering rate below T^ due to the opening of the supercon- 
ducting gap, but it is hard a priori to guess why some- 
thing similar should happen in the presence of an elastic 
scattering rate of order lOOmeV. The second problem is 
that recently increased momentum resolution^ and the 
use of different photon energiesiSiiiiiSii^ has resolved a 
bilayer splitting which has its maximum effect near the 
(tt, 0) point. Some of the previously observed "elastic 
broadening" is therefore certainly due to this as well as 
to pseudogap effectsi^, but exactly how much is not clear. 

On the other hand, one cannot ignore the out-of-plane 
disorder. The BSCCO material is thought to be doped 
by excess oxygen in the SrO and BiO planes, and even 
the best single crystals are believed to contain signifi- 
cant amounts of cation switching and other out-of-plane 
defectsi^. It is therefore reasonable to assume that quasi- 
particles moving in the Cu02 planes of this material must 
experience a smooth potential landscape due to these 
defects, and useful to pursue the question of the effect 
of this type of scattering in the superconducting state. 
In fact, fitai^ to Fourier transformed-scanning tunnelling 
spectroscopy measurementsii on similar samples to those 
used in the ARPES studies have recently been shown to 
require both a strong (near-unitary limit) scattering com- 
ponent, attributed to native defects in the Cu02 planes, 
as well as a weaker, smooth scattering potential com- 
ponent attributed to defects away from the plane. Re- 
cently, Markiewicz has also attempted to relate STM and 
ARPES data in the superconducting state assuming a 
smooth potentialii^ 

In this work we model the complex collection of out- 
of-plane defects with a simple set of impurity potentials 
with finite range, and find a number of surprising re- 



suits. The first is that the elastic scattering rate indeed 
"collapses" in the superconducting state in the strong 
forward scattering limit, leading to a sharp spectral fea- 
ture everywhere on the Fermi surface except at the nodal 
point itself. The second is that as one goes away from 
the Fermi surface, this feature disperses away as expected 
for a quasiparticlc peak only if the system is sufficiently 
clean; otherwise it remains pinned to the gap edge. This 
would appear to explain the apparent lack of quasiparti- 
cle dispersion in older samples. In some situations both 
a nondispersive gap edge peak and a dispersing quasi- 
particlc feature can be simultaneously observed. While 
we do not attempt a direct fit to experiment, we show 
that assuming a rather simple and physically motivated 
model for the one electron self-energy, which combines 
the above description of forward clastic scattering with 
strong pointlike elastic scattering and spin fluctuation in- 
elastic scattering, allows us to calculate a spectral func- 
tion which appears to reproduce many of the qualitative 
features of current ARPES data on optimally to over- 
doped cuprates. It is our hope that the ideas presented 
here can help unravel some of the mysteries surrounding 
the behavior of what is generally called the supercon- 
ducting "quasiparticlc peak" , and allow a more accurate 
description of the actual propagating excitation. They 
will also have important immediate implications for other 
bulk properties of the superconducting state. 

The paper is organized as follows. Section II discusses 
general effects of extended impurities scattering on the 
ARPES spectrum in both the normal and superconduct- 
ing states. In Section III, we examine the effect of adding 
an isotropic elastic scattering due to unitary scatterers as 
well as the inelastic scattering due to spin fluctuations. 
Finally, section IV contains our comparison with existing 
data, conclusions and plans for future work. 



II. ELASTIC SCATTERING 

A. Normal State 

We first consider a model system where scattering oc- 
curs only because of disorder. For simplicity, we will 
assume that the most important feature of the poten- 
tial due to out-of-plane impurities experienced by elec- 
trons moving in the Cu02 planes is its finite range k~^. 
We therefore model a single impurity simply as a term 



Vir) = Voe- 



V, 



2ttkVq 



kk' 



((k-k')2+K2)3/2' 



(1) 



where Vq sets the strength of the potential and k"^ is its 
range. Note that k and k' are only defined up to a re- 
ciprocal lattice vector. The self-energy due to many such 
impurities gives rise to an elastic broadening of quasi- 
particlc states which depends upon the position of k in 
the Brillouin zone. For weak scattering, the Born ap- 
proximation for the retarded self-energy associated with 
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FIG. 1: The Fermi surface corresponding to the band efc given 
by Eq. Q with f'/f — —0.35 and fi — —1. The energy distri- 
bution curves of A(k,uj) will be discussed for k — kjv (nodal), 
k = kA (antinodal), and various k values along the (0,7r) to 
kA cut shown by the solid points. 



a random distribution of n/ impurities per unit area has 
the usual form in the normal state 



E(k,o;)=n/V|VukfG°(k',u.), 



k 



(2) 



where oj is understood to include an infinitesimal positive 
imaginary part, and the superscript on the single-particle 
Green's function G^ = (w — ek + aO^^ indicates that at 
first we ignore self-consistency, i.e. G*^ does not depend 
on E. 

To model the electronic structure we take a simple 
near-neighbor hopping t and next-near-neighbor hopping 
t' , such that 



Ek — — 2t(cos fcj. 4- cos ky 



At' cos fcj. cos ky ~ fi 



(3) 



and set t'/t = —0.35 and ji/t = —1. Note that with 
this choice of parameters there is a van Hove singularity 
with a peak in the total density of states at — 0.4i. The 
Fermi surface for these parameters is shown in Fig. ^ it 
is similar but not identical to the Fermi surfaces found 
for both YBCO-123 and BSCCO-2212 by ARPES. 

As the range of the potential k~^ increases, the scat- 
tering of a quasiparticlc from k to k' becomes peaked 
in the forward direction. As shown in Figure El when k 
is close to k' and both are not too far from the Fermi 
surface, we may parameterize them as 



k = k^ + kj_kj_ 
k' = ki? + qii + k'j_k'j_ 



(4) 

(5) 



where q = k — k' is the momentum transfer and q|| its 
component parallel to the Fermi surface. The unit vec- 
tors k± and k'j^ are the projections of k and k' onto 
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FIG. 2: Geometry for the forward scattering process in which 
a quasiparticle scatters from k to k'. 



the Fermi surface, respectively, such that, e.g., Ck' = 
vp{'k')k'j_. The imaginary part of the retarded self-energy 
(O becomes 
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FIG. 3: The elastic scattering width r(k) = -ImE(kF,0) 
in units of t plotted around the Fermi surface with 9 = 



X cLlT I rCy / Kix 



Results are shown in Fig. 2a-d for 



10, 5, 2, and 0.5 respectively with n/|Vbp adjusted so that 
Ymax = — Im E(kA,0) is equal to 0.2f. The dashed curve is 
proportional to Vp^{9) and one sees that as k decreases and 
the scattering peaks in the forward direction, r(k) varies as 

v-/{e). 
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37rn,Vo^ 
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Eq. ((HJ shows explicitly that in the limit of small k, 
the self-energy becomes more and more sharply peaked 
"on the mass shell" lo = Ck- This is a generic feature of 
long-range potentials. For example, on the Fermi surface 
k = \ip at cij = 0, 



-j:"ikF,o) = ro(,kF) 



37rn,Vn^ 

\vF{'kF)\K^ 



(9) 



In Figure we show how this angular dependence oc 
I/Ii'IcfI is approached by the exact result Q as k de- 
creases and the range of the scattering is increased. In the 
figure and in what follows in this section, we will use the 
result with n/|VbP chosen such that ro(kyi) = 0.2i 
corresponding to ro(kyi) ~ 30 meV for t = .15 eV. 
From Fig. 3 we see that in the forward scattering limit, 
ro(kjv) ^ ro(kA)/1.4 ~ 0.14i. We stress that the pre- 
cise dependence of the forward elastic part of the self- 
energy © on K and on momentum depend on the details 



of the Fermi surface shape and impurity scattering po- 
tential. However, we do not expect qualitative features 
of the resulting spectra to be affected. We note further 
that the absolute magnitudes of the parameters chosen 
arc roughly consistent with the 0(10%) weak scatterers 
of strength V^ of 0{t) and range of O (l-2a) extracted 
from FT-STS data in Refpi&. In Section III, we will show 
results for various values of ro(kyi) and k. 

Thus far we have not considered the effect of self- 
consistency, i.e. replacing G" in Eq. Q by G(k, lo) = 
[(G°)~^ — S]"^. In the normal metal one is used to 
ignoring this distinction, as the self-consistent solution 
for pointlike scatterers can be shown to be identical to 
the non-self-consistent one up to corrections of order 
{oj/Ef)'^, where Ef is the Fermi energy. If the scat- 
terers have a finite range, however, this argument breaks 
down and self-consistency becomes important. As seen 
in Figure ^ the effect of self-consistency is to reduce 
the frequency dependence of S"(k, w) induced by elec- 
tronic structure; in particular, the van Hove singularity 
at uj = — 0.4t is eliminated. This may account for the 
complete absence of van Hove spectral features in STM 
and other tunnelling experiments on BSCCO. We note 
further that the scattering rate is cut off in the forward 
scattering case k = 0.5 when \lu\ > vfk, reflecting the 
fact that if the electron's momentum can only be shifted 
by a small amount ~ k, the allowed energy transfer is 
also restricted. 
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FIG. 4: Scattering rate r(k, tj) = — ImE(k, tj) for k on the 
Fermi surface at points A (left) and N (riglit), for two differ- 
ent values of forward scattering parameter k = 5 (top) and 
0.5 (bottom). Here, as in Fig. El ntVo has been chosen to 
give r(kA) ~ 0.2i. Dashed curves: non-self consistent Born 
approximation; solid curves: self-consistent Born approxima- 
tion. 



B. Superconducting State 

1. Elastic Self-energy 

To describe the superconducting state, we assume 
a BCS d-wave order parameter corresponding to pure 
nearest-neighbor pairing. 



Ai 



— (coak^ 



cos ky) 



with Ao/i = 0.2. The full matrix Green's function in the 
presence of scattering in the superconducting state is 



Gik,uj) 



UJTg + £kT3 + AkTl 



(10) 



where tD = w — Sq, Ck = Ck + S3, Ak = Ak + Si, and the 
Sq, are the components of the self-energy proportional to 
the Pauli matrices r^ in particle-hole space. If we first 
assume the simplest case, that the scattering is entirely 
elastic and weak, we may approximate the self-energy in 
the Born approximation similar to ^ as 



ni 



Z^ 



|^kk'|'T3G°(k',C^)T3, 



with Nambu components 

Eo(k,c^) = niY,\V{k,k')\^ 

k' 

I]3(k,t.) = n, ^|y(k,k')P 

k' 

and 



C.2 - el - ^l, 



Ck' 



A^ 

^k' 



(11) 
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FIG. 5: The self-energy terms -Im Eo(k, a;), Im Ei(k, oj), and 
-Im E3 (k, ijj) in the superconducting state at T = for k = k^ 
(top) and kiv (bottom), for «: = 5 and 0.5 and the same band 
and scattering parameters as previously used. Here Afe = 
Ao (cosa; — cos ky)/2 with Ao = 0.2f . 
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We will calculate the self-energies both non-sclf- 
consistently, as in Eq. (|ll|l . and self-consistently by re- 
quiring that S[G°] -^ S[G]. In Figure [S] we show the 
variation of these self-energy components with energy at 
k = kjv and k = k^. Again the van Hove singularity 
is washed out in the self-consistent evaluation, and it is 
furthermore noteworthy that the E3 component becomes 
quite small in the forward scattering limit. While in Fig- 
ure |S1 gaps appear in the Sq near the A point for small 
K, the values of Eg and Si near the gap edge to > Ak 
are large, of order several times the hopping t. Quasipar- 
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FIG. 6: Scattering rate V^i{V.,ijj) vs. u) for k = kA(left) and 
kAT (right) in the superconducting state at T = 0, for k = 5 
(top) and K = 0.5 (bottom). Here ro(kA) = 0.2i. 



ticle properties near the Fermi surface are determined, 
however, by particular combinations of the Nambu self- 
energy components. As can immediately be seen from 
the denominator of (|1U|I , the total elastic scattering rate 
broadening the quasiparticle state of energy u will be 



Te^k) ~ -Im So(k,c^) + — Si(k,c^; 



(15) 



provided one can neglect S3 (see Fig. EJ. In Figure El 
we see that as uj goes to --A(k^) for k = k^, Y^i is sup- 
pressed by the near cancellation of the two components 
in H15|) when k becomes small. To obtain some insight 
into the physical origins of this cancellation, we derive 
approximate analytical forms for k <C 1 following the 
discussion of the normal state above, leading to 



Eo(k,w) -- -ro(k_F) 

E;'(k,c^) ^ ro(ki.) 

E^'(k,c.) ~ 0, 



v4J^^A 
AkSgntj 



(17) 

(18) 
(19) 



but are strongly supressed due to energy conservation 
when \\/ijS^ ~ A^ — ek| becomes greater than kvf-, as one 
may observe in Figure [S] These equations are now iden- 
tical in form to those expected for an s-wave supercon- 
ductor (even when self-consistency is included). All pecu- 
liarities of the d-wave state which result from momentum 
averaging over the Fermi surface have disappeared. We 
therefore expect a 'priori to recover Anderson's theorem, 
the insensitivity of bulk thermodynamic properties to 
nonmagnetic scattering2£. The physical reason for this is 
clear: as seen in Figure|21 for small k and k a distance at 
least K from the node, small angle scattering cannot mix 
order parameters of different signs, and therefore cannot 
break Cooper pairs. The analogy with the s-wave su- 
perconductor and relation to Anderson's theorem in the 
context of ARPES is discussed further in Appendix 1. In 
addition, the slowing of the rate of T^ suppression due 
to disorder as scattering becomes more anisotropic was 
treated some time ago by several authors^. These works 
were motivated by T,^ suppression rates in the cuprates 
which appear to be 2-3 times slower than predicted by 
the classic Abrikosov-Gorkov formula^i appropriate for 
for pointlike isotropic scatterers. The suppression of T^. 
near the pure forward scattering limit has been discussed 
recently in detail by H.-Y. Kec»2,. 

In the forward scattering limit where H17I19|I hold, the 
effective elastic scattering rate (|15|) becomes 



re/(k,tj) ~ ro(kF^ 
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A^ for the 



Here Sq — |aj|, — AkSgntj, and v^'gnu^fu 
Nambu components a = 0, 1 and 3 respectively, and 
ro(k) is given by Eq. Q. Note that E3 vanishes on 
the Fermi surface ek = in this limit. In order to gain 
some further intuition for these expressions, we special- 
ize to the case where the momentum k is close to the 
Fermi surface and the energies oj are small, such that 
\(e.k ± •\/w^~-"Aj)/KWi?(k)| ^ 1. The self-energies may 



For k along the nodal direction, the elastic broadening 
in the superconducting d-wave state is equal to its value 
in the normal state. However, for k in the antinodal re- 
gion, the broadening vanishes as w ^ Ak and approaches 
the nornral state value only when lo becomes large com- 
pared with Ak- The elastic contribution re/(k^,a;) to 
the broadening at the antinodal \i.A point versus lo is 
shown in Fig. |S1 Physically, the individual contribu- 
tions to the normal I]o(k, lj) and anomalous Si(k, w) 
self-energies are both enhanced by the density of states 
factor {up' — A^(k))~2 (Fig. EJ. However, the normal 
contribution describing the scattering out of state k into 
k' is compensated by the anomalous contribution scat- 
tering into k from the pair condensate. This gives rise 
to the suppression of the elastic scattering seen in Fig. 
relative to Fig. [Sjas to approaches — |Ak^| from below. 



2. Spectral Function 

The near-cancellation of the two Nambu components 
of the self-energy near the gap edge in the forward scat- 
tering limit leads to a dramatically reduced elastic scat- 
tering rate in the superconducting state, which sharpens 
the spectral features of quasiparticles which are not too 
close to the nodes. Within the current model where we 
continue to neglect inelastic scattering due to electron- 
electron interactions, we now turn to the 1-electron spec- 
tral function measured by ARPES. In the forward scat- 
tering limit, with the self-energy given by Eqs. (fT7jl - (f11^ . 
one obtains a result for the Green's function previously 
discussed by Markicwicai^, 

UI^K, CJj_ -— 5 . [^IL) 

(cj2- A^)z(k,w)^-e^ 



Here z(k,w) = 1 + iro(k)sgnti;/-\/w2 — A^.. The electron 
component of the spectral function is then 



yl(k,w) = ImGii(k,w) 

TT 

1 tjz(k,w)-fek 

= Im 
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It is useful to consider a few special cases of (|22f) more 
closely. In particular, on the Fermi surface ek = one 
has the simple expression 



^(kF,w) 



ro(k) 



w 



while near the gap edge, e.g. lu < — |Ak|, 



:, (23) 



TT el + r2(k) ^IJ?3a 



(24) 



At the nodal point k^v, where the gap vanishes, the 
spectral weight is given by the simple Lorentzian form 



A{\<iN,uj) 



ro(kAf)/7r 

'2+ro(kA.)2' 



(25) 



and at low temperatures where the elastic scattering is 
dominant one can determine ro(kjv). However, for k at 
the antinodal point k^ such that — Ak^ — doj < to < 



-Ak. 



A{]f.A,Oj)'^ 



1 A 



kA 



'^ ro(kA) ./^2_A2^ 



(26) 



where the "width" 5lj depends upon the ratio 
AkA/ro(kA)- If Ak^ is large compared with ro(k^), the 
width 5uj ~ ro(kyi)2/(2Ak4). If one integrates ^(k^,^) 
from — Aki — 5lo to — Ak^ to define a "peak intensity" 
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FIG. 7: Comparison of the T = self-consistent normal state 
(dashed) and superconducting (solid) one-electron spectral 
functions for k = 5 (top) and 0.5 (bottom) at k = k^ (left) 
and kiv (right), with ro(kA) = 0.2f. 



which is independent of Ak^/Fo(k^). However, when 
the system is sufficiently dirty such that Ak_4 < Fo(k^), 
the falloff of A{\i.A, io) as lo decreases below — Ak_4 varies 
as (w2 — A^^)~^/2^ jn ti^js case, the scale is set by Ak_4 
and if one takes 5uj ~ Ak4, the peak intensity varies 
as /(k^) ~ Ak_4/Fo(kA). This is quite different from 
the usual BCS quasiparticlc result which is proportional 
to the quasiparticlc renormalization factor 2;(k^) times 
a coherence factor which is 1/2 on the fermi surface. It 
should be possible to test the foward scattering scenerio 
by comparing the variation of /(k^) with Ak^/Fo(k^). 
Numerical results for the normal and superconducting 
spectral weights for w < using the self-consistent ver- 
sion of the self-energy H12I14|I with tj ^ w, etc., for the 
model impurity potential (1) arc shown in Fig.Qfor k = 5 
and 0.5. For ti = 0.5, the scattering is predominantly 
forward and in the superconducting state at T = 0, the 
spectral weight for k^ is seen to sharpen at — Ak^ de- 
spite the fact that the normal state broadening F(kyi) 
is of order the peak position. For k = kjv, the nodal 
spectral weight has the expected Lorentzian form with 
a width set by r(kAr). For k = 5, corresponding to a 
more isotropic scattering, the spectral weight for k = k^i 
broadens. In the limit of isotropic impurity scattering, 
the El component of the impurity self energy vanishes 
and y4(kyi,a;) has a Lorentzian-like broadened peak. 



III. ISOTROPIC ELASTIC AND INELASTIC 
ELECTRON-ELECTRON SCATTERING 

In addition to the forward scattering by out-of-plane 
impurities, in this section we consider unitary limit 
isotropic clastic scatterers as well as inelastic electron- 
electron collisions. In fact, it is the momentum and frc- 




FIG. 8: Schematic depiction of the various contributions to 
the scattering rate in tire superconducting state at T = 0. 
First panel; elastic forward scattering for k = Ica and kjv; 
Second panel: isotropic unitarity limit impurity scattering 
due to ~ 0.2% impurities; Third panel: spin fluctuation in- 
elastic scattering rate interpolated from Ref^ 



quency dependence of the latter interaction about which 
one hopes to learn more from the ARPES spectrum. Here 
we proceed phcnomcnologically by writing 



^tot 



^elj + ScLu + Sincl- 



(28) 



The first term is the contribution from elastic (quasi-) for- 
ward scatterers we have discussed in Section II. The sec- 
ond term represents the effect of unitary scatterers (possi- 
bly Cu vacancies) with concentration roughly n„ ^ 0.2% 
observed as zero-bias resonances in STM experiments. It 
will be treated as usual in the self-consistent T-matrix 
approximation, 



5eL- 



EkG(k, 



:^o- 



(29) 



In the normal state we find a scattering rate of r„ ~ 
10^'^t, leading to an impurity bandwidth ^u ~ V^u^o ^ 
10~^t (of order 1 to 2 meV). We note that the width of 
the resonance observed in STM is in fact roughly 3 times 
this number, consistent with the self-consistent T-matrix 
calculation^; impurity resonances are visible up to about 
10 meV in experiments. Nevertheless the isotropic part of 
the elastic scattering appears to have a relatively insignif- 
icant effect on the ARPES spectral function, as shown 
below. 

For the imaginary parts of the inelastic self-energies, 
we will use numerical results obtained from a spin- 
fiuctuation calculation of the quasiparticle scattering. 
Following the argument leading to (|15|l . we define an ef- 
fective inelastic scattering rate 

\k,iu) + ^j:r\h,u)). (30) 



Find (k, uj) = -Im S 



In spin-fluctuation calculations of Find (k, uj) , it was 
foun d"^«- ' n that at the nodal point, at low temperatures, 
the scattering rate initially increased as the third power 
of Lu or T depending upon which is larger. At other k 
points on the Fermi surface, the scattering rate varies ap- 
proximately as the third power of this energy measured 
relative to A^. The reduction of the inelastic scattering 
rate at low excitation energies reflects the suppression of 




t« 



e. 4 



-0.25-0.2-0.15-0.1-0.05 0.05 

(0 

K=0.5 



I ,-11 



-0.25-0.2-0.15-0.1-0.05 0.05 

CO 



FIG. 9: Finite temperature spectral function at the antinodal 
point A on the Fermi surface multiplied by the fermi func- 
tion, A{\i.A,i^)f{i^) vs. UJ including the full model self-energy 
as described in Section III. Results for k = 2 and 0.5 with 
ro(k.4) = 0.2t are shown. 



the low energy spin-fluctuations due to the opening of the 
d-wave gap. Here, we will use a numerical interpolation 
of the w-and T-dependent Find obtained from Refi2^. 



A. The Antinodal Spectrum 

The various contributions to the self-energy are 
sketched in the schematic diagram shown in Fig. |H1 We 
will use the parameters discussed above to set the mag- 
nitudes of Eind and Sd,u. Then we will consider various 
values of k and Fo(kA), characterizing the forward elastic 
scattering. While the effect of Sind is already contained 
in the intrinsic Tc and the suppression of T^ due to Sci,u 
is negligible for the parameters considered, this is not 
necessarily the case for the forward elastic scattering. As 
discussed by Kee^, the suppression of Tc varies as k^, so 
that for small values of k such as k = 0.5, the suppression 
of Tc is negligible for the scattering rates Fo(kj4) which 
we will consider. However, at larger values of k this is not 
the case, so that the Tc shift or the lack of shift implies 
constraints on k and Fq. For the moment, however, we 
will ignore this and simply compare k = 0.5 with k = 2 
for three different scattering rates Fo(kyi) = 0.2, 0.1, and 
0.05 in units of t corresponding to Fo(kA)/Ao = 1, 0.5, 
and .025. 

In Fig.|31we show the temperature dependence of the 
energy distribution curves A(k, w) /(w) for k at the antin- 
odal point k^. Here /(w) is the Fermi function. In the 





FIG. 10: A(k,uj) vs. u! for k = 2 and 0.5. Results are given 
for the k points as indicated along the (0, n)- A cut in Fig.Q 
The disorder levels correspond to r()(kA)/Ao = 1, 0.5, and 
0.025. Note the spectra for different k points have been offset 
for clarity. 



superconducting state there is a square-root-like behav- 
ior as uj approaches — |Ak^| from below. This should be 
contrasted with the broad Lorentizian peak in the nor- 
mal state which is cut off by /(w). While one could 
understand that inelastic broadening would diminish as 
the temperature is lowered, the asymmetric, one-sided, 
square-root-like sharpening of the spectrum in the su- 
perconducting state is a consequence of forward elastic 
scattering as discussed in Section II. The peak intensity 
at k^ should scale as Ao(T). 



B. Quasiparticle Dispersion Near the Antinodal 
Point 



In a clean superconductor, there is a peak in the spec- 
tral function A(k, w) at the quasiparticle pole lj = E]^ = 
•\/e^ + A^. In particular, as the momentum moves along 



the cut from (kyi,7r) to (0, tt) shown in Fig. 1, one ex- 
pects to see a dispersion of this peak to higher ener- 
gies. However, if the forward elastic scattering strength 
ro(kA) ^ Ak_4 , then the peak in A(k, oj) remains at — Ak 
rather than dispersing. Figure [TUl shows plots of A(k, a;) 
for different values of k between the M and A points for 
K = 2 and 0.5, and several values of the scattering rate 

ro(kA). 

As samples improve, there is a natural tendency in this 
model for the spectrum for k not too far from the antin- 
ode to cross over from one characterized by a nondisper- 
sive peak at Ak in the dirty limit where Fq ~ Aq to 
one characterized by a dispersive quasiparticle peak at 
i^k when Fq is small compared to Aq. This crossover 
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FIG. 11: Finite temperature spectral functions at the nodal 
point N on the Fermi surface, A(kjv,aj)/(aj) vs. ui including 
the full model self-energy as described in Section III for k = 2 
and 0.5, with Fofkiv) = 0.14t. 



is due to the way in which the forward elastic scatter- 
ing rate for a d-wave superconductor is reduced as the 
gap edge is approached and is analogous to the same ef- 
fect discussed analytically in Appendix 1 for an s-wave 
superconductor. In a system with Fq ^ Aq, no quasipar- 
ticle peak is observed, but a sharp feature does appear at 
— Ak, representing simply the spectral weight in the over- 
damped quasiparticle peak piling up at the gap edge as in 
the s-wave case. Only when Fq becomes small compared 
to Aq does one see a true quasiparticle peak dispersing 
as —Ek- In the most strongly forward scattering case, 
K = 0.5, one can see that, depending on the strength of 
the scattering rate, one can have simultaneously a broad- 
ened dispersing feature as well as a gap edge feature. It 
is tempting to speculate that this phenomenon is related 
to the peak-dip-hump features observed generically be- 
low Tc in cuprate ARPES experiments, but we have not 
yet explored this issue in detail. 



C. The Nodal Spectrum 

In Figure [TTl the corresponding EDO's for the nodal 
k^ point are shown. In this case, the T — spec- 



tral function is a Lorcntzian centered at the Fermi level, 
whose width is limited essentially by the elastic scatter- 
ing. For the band structure parameters we have chosen, 
ro(kAr) ^ ro(kA)/1.4, so that ro(kAr) = Q.Ut. Results 
in Fig.^]are shown for several different values of ro(kjv). 
At finite temperatures the peak is further broadened by 
inelastic processes. Defining a width of the asymmetric 
EDC's is quite difficult, a natural result of the strongly 
w-dependent self energy in the current approximation. 
In fact, even the shape of the EDC curve is difficult to 
compare directly to experiment, since it can be qualita- 
tively changed by a small amount of averaging along the 
Fermi surface due to the angular resolution of the detec- 
tor. If one goes a short angular distance away from the 
node along the Fermi surface, one finds at low T not a 
smooth Lorcntzian for uj < 0, but rather a square root 
singularity with a small local gap Ak. Averaging over a 
small k region will therefore make the "leading edge" of 
this spectrum appear much sharper. 

The intrinsic T-dependent broadening is therefore 
most clearly seen in the momentum distribution curves 
(MDC's), which in this work (see Eq. H25|) ) arc simply 
Lorentzians centered at the Fermi surface. In Figure [T^ 
we show that, within this model, the rate of sharpening 
of the nodal MDC's indeed increases somewhat in the 
superconducting state as the temperature is lowered, but 
that the widths saturate at a value determined by the 
elastic scattering. The temperature dependence of Afc 
shown in Fig. E|is in disagreement with the results of 
Valla et ali^, who reported a linear MDC width depen- 
dence on temperature even in the superconducting state. 
However, the type of lo and T dependence wc find would 
appear to be consistent with nodal EDC linewidth mea- 
surements of Kaminski et al^. 
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IV. CONCLUSIONS 

The main point of this work is to raise the possibility 
that some of the ARPES observations on optimally or 
overdoped samples which appeared to be in conflict with 
the BCS results may in fact fit within the BCS frame- 
work when the effects of forward elastic scattering are 
taken into account. Thus, while forward elastic scatter- 
ing can be responsible for some of the anomalous width of 
the spectral function measured above Tc, this need not 
be in conflict with the observation of a sharp spectral 
feature in the antinodal region below Tc- Furthermore, 
if ro(k^) is larger than A(k^), the intensity associated 
with the area under this feature varies as A(kyi)/r(kyi), 
which depends upon the doping and temperature. At 
the same time, the spectrum at the nodal point can ex- 
hibit a Lorcntzian behavior with a width that evolves 
smoothly through Tc and then partially narrows as the 
inelastic scattering is suppressed by the opening of the 
gap. We have seen that these effects occur because of 
a phenomenon similar to Anderson's theorem which ap- 
plies for much of the Fermi surface of a d-wave supcr- 



FIG. 12: a) Finite temperature spectral function near the 
nodal point N along a cut perpendicular to the Fermi surface 
at N, A(k,Li} = 0) vs. k, including full model self-energy with 
K = 0.5 and ro(kiv) = O.Ut. b) The half-width oiA(k,uj = 0) 
at half-maximum, Afc, plotted vs. T/T^. 



conductor in the forward scattering limit: if scattering 
is sufficiently peaked in the forward direction, the scat- 
tering process does not mix states with different signs of 
the order parameter, and no pairbreaking occurs. Suf- 
ficiently near the nodal point, however, the sign change 
always takes place, implying that the full width of the 
quasiparticle due to elastic scattering is recovered at the 
nodal k^r. 

We have discussed simple approximate forms as well as 
fully self-consistent numerical calculations for the elastic 
self-energies. In the overall model for the self-energy, 
we also included an approximate treatment of electron- 
electron collisions, as well as unitary "native defect" scat- 
terers observed to be present in the BSCCO-2212 mate- 
rial in STM experiments. The energy distribution curves 
found in the self-consistent calculations using this model 
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show no influence of the van Hove singularity, and dis- 
play an asymmetric shape with a rounded square root 
peak near the local gap edge. Widths and temperature 
dependences can be obtained which appear comparable 
to experiment, but it is difficult to compare directly be- 
cause of the unknown ARPES background signal and bi- 
layer splitting, which was not included here. We also 
discussed the temperature dependence of the momentum 
distribution curves in the superconducting state at the 
nodal point. It was found to have a Lorentzian shape 
which narrows with decreasing temperature in the su- 
perconducting state, due to the suppression of inelastic 
scattering, saturating at a value determined by the elastic 
scattering. 

Along the momentum cut (kyi, tt) to (0, tt), if ro(kA) is 
small compared to Ak, we find a dispersing quasiparticlc 
peak at a; — — £'fc. However, when ro(kyi) is comparable 
with Ak, the maximum response occurs for uj ~ — Ak, 
does not disperse as k moves away from the Fermi level, 
and can still represent a sharp spectral feature. This is 
similar to some of the older data on the BSCCO-2212 
system, which has not been explained at this writing and 
might have been taken for some new type of dispersion- 
less excitation in the superconducting state. However, as 
we have shown, it is simply the consequence of forward 
elastic scattering which is suppressed at the gap edge in 



the superconducting state. 

Here we have focussed primarily on the effects of for- 
ward elastic scattering on the BSCCO-2212 ARPES spec- 
trum, and found that a model of several per cent weak 
out-of-plane scatterers with a range of order one lattice 
spacing, similar to current models of Fourier transform 
STM measurements, can explain many qualitative fea- 
tures of the ARPES data. We believe that the unique 
way in which this material is doped and disordered using 
current crystal growth techniques endows it with an effec- 
tive disorder potential which strongly influences the low- 
energy quasiparticlc properties in a way which is charac- 
teristically different from the much cleaner YBCO sys- 
tem, for example. This picture should have important 
and calculable consequences for other superconducting 
properties, such as microwave conductivity, which wc ex- 
plore elsewhere. 
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APPENDIX A: ONE-ELECTRON SPECTRAL 

FUNCTION IN A DISORDERED S-WAVE 

SUPERCONDUCTOR 

Insight into the spectral function of a rf-wavc supercon- 
ductor with k near the antinodal point and reasonably 
forward scattering can be gained by examining the spec- 
tral function of an s-wave superconductor with disorder. 
Although quite simple, we are not aware that the spectral 
weight for an s-wave superconductor with isotropic im- 
purity scattering has been discussed elsewhere. For sim- 
plicity, we restrict our consideration to isotropic, weak, 
nonmagnetic scatterers. The self-energies in the Born 
limit for a system with particle-hole symmetry and an 
isotropic order parameter A are 



Eo - n,1/2^- 



k ^'-el-A^ 



-zriS^ (Al) 



El = n,V^J2- 



A 



= iT- 



sgnw' 



uj"^ -el- A2 Vu^- 1 



E, = 0, 



(A2) 
(A3) 



where u = Co/ IS. and u' is Re u. The self-consistency 
equation is then 



cj — Eo (tj/A)\A?--T + Jsgnu'Tw /. .^ 
A + Ei ^ Vu^-l + isgnuT ' 

with r = r/A, which has the solution 



In Figure [TSI we plot the spectral function ^(k, cj) = 
-Im G^ 1 fk. Lo) liT for u; < for various values if k on and 
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FIG. 13: j4(k,(i)) for an s-wave superconductor vs. c^j/A for 
r/A = 0.1 for 3 values (solid, dashed, dotted) of tk/A = 
0,-1,-2. 



near the fermi surface. Since iT j \/uP' — A^ is pure real 
for Icjj/A < 1, the spectral weight vanishes for energies 
|a;| < A the (unrcnormalized) gap. However, as lo ap- 
proaches —A from below, the spectral weight for k = kp 
diverges as 



« = -^ = — . A5 

A A ^ ' 

This means, as is well known, that the density of states 



p{uj) = Re 



CO sgn uj 
VCu^ - A- 



Vuj^ - A2 • 



(A6) 



is not changed by this type of impurity scattering and 
there is no rcnormalization of the momentum-integrated 
thermodynamic properties. This is the essence of "An- 
derson's theorem"—. On the other hand, one might ex- 
pect that the spectral weight A(k, uj) for a quasiparticle 
of momentum k should be broadened by disorder since 
the scattering mixes different momentum states. This is 
the case in the normal state where the spectral weight in 
the presence of impurity scattering becomes a Lorentzian 
of width 2r. To determine what happens in the supercon- 
ducting state we need the 11 component of the Nambu 
Green's function 



Gii(k,c^) 



1 



iFsgn UJ 

7z 



-A2 



fifc 



iFsgn UJ 



(A7) 



(^2_A2)-e| 



A(kF,uj) 



A 



ttF V^ 



(A8) 



As seen in Fig. 1131 this square root singularity gives way 
to a more symmetrically shaped dispersing quasiparticle 
peak as k moves away from k^?. However, for k not too 
far from the Fermi level, a residual square root singu- 
larity at w = — A remains. In Figure 1141 which shows 
the dependence of the spectra on disorder, we see that 
away from the Fermi level the strength of the structure 
at w = — A increases with disorder. In fact one can eas- 
ily show that when |efc| > F, the spectral weight as uj 
approaches —A varies as A{'k,uj) ~ ^-% , / ^ . Thus 

in cleaner systems, the anomalous peak at a; = —A dis- 
appears and one has just the expected quasiparticle peak 

at a; = —Ek- 
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FIG. 14: A(k,u}) for an s-wave superconductor vs. uj/A for 
Ek/A = -1 for 3 values (solid, dashed, dotted) of T/A = 0.01, 
0.2, 0.5. 



